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We study charmed baryon resonances which are generated dynamically within a unitary meson- 
£SJ ' baryon coupled channel model that treats the heavy pseudoscalar and vector mesons on equal 

footing as required by heavy-quark symmetry. It is an extension of recent SU(4) models with t- 
channel vector meson exchanges to an SU(8) spin-flavor scheme, but differs considerably from the 
SU(4) approach in how the strong breaking of the flavor symmetry is implemented. Some of our 
dynamically generated states can be readily assigned to recently observed baryon resonances, while 
others do not have a straightforward identification and require the compilation of more data as 
well as an extension of the model to d-wave meson-baryon interactions and p-wave coupling in 
the neglected s- and ii-channel diagrams. Of several novelties, we find that the A c (2595), which 
emerged as a ND quasi-bound state within the SU(4) approaches, becomes predominantly a ND* 
O l' quasi-bound state in the present SU(8) scheme. 
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I. INTRODUCTION 

> 

^v ' Iii hadron physics establishing whether a resonance has the usual qq or qqq structure or better qualifies as being 

\^ , dynamically generated through the multiple scattering of their meson and/or baryon components, has kept an active 
0^ ' topic of research. In recent years, the introduction of unitarity constraints in coupled channels methods has led to 
extensions of the chiral effective theories that can be applied at much higher energies and in the vicinity of resonances 
[l|. In the particular case of baryons, several resonances are firm candidates to be identified with states generated 
dynamically from the interaction of mesons of the 0~ octet Goldstone bosons with the l/2 + octet ground state 
baryons @, |, i, i, 0, 0, H, i, 0, M MJIlJliJEJlE EDI El GJ, EH ■ Then interesting predictions such as the 



two-pole nature of the A(1405) [3, [lj, [2l|, [22, [23, uA, l25| have found experimental confirmation [2fj, [27|, as discussed 
in Ref. j28| . Later a number of works has also been devoted to the study of dynamically generated J p = 3/2~ states 
[29L [30l. T3ll. [321 . 133 . 134 . I351 . [3a . [37I [381 [39j |. Early works considered the interaction of pseudoscalar 0~ mesons with 
the baryons of the 3/2+ decuplet. The incorporation of vector mesons into the above-mentioned coupled-channel 
picture was pursued in studies of axial 1 + meson resonances of Refs. [40, |41|. b ut until very recently those degrees 
of freedom have not been considered in the baryon- meson sector. In Ref. [37j | the (83,1a) 1 nonet of lowest-lying 
vector meson (p(770), if (890), -R"(890),u(780) and </>(1020)) degrees of freedom was incorporated into the Weinberg- 
Tomozawa (WT) meson-baryon chiral Lagrangian by using a scheme which relies on the SU(6) spin-flavor symmetry. 
The corresponding Bethe-Salpeter (BS) equation successfully reproduces the previous SU(3)-flavor WT results for 
the lowest-lying s- and d- wave, negative parity baryon resonances (iV(1535), iV(1650), A(1390), A(1405), A(1520), 
A(1690), A(1670), S(1620), S(1690), • • • )■ It also provides some information on the dynamics of the heavier ones such 
as A(1800) or A(2325) which have sizable NK* and Aw couplings, respectively |38l.l39l|. 

Recently the attention has turned towards studying resonances with charm degrees of freedom, motivated by the 
discovery of quite a few new states, as reported by the CLEO, Belle and BaBar Collaborations [42. lia. H. lia. I461. |47| . 
At first sight, the physics in the charm C — \ sector bears a strong resemblance to the phenomenology seen in the 



1 Here and below we use the notation /i2j+i, for the SU(3) multiplet fi of spin J. 



KN dynamics upon replacing an s quark by a c quark 2 . This is reinforced by an apparent similarity between the 
two I — s-wave A(1405) and A c (2595) resonances. The mass of the former lies in between the 7rS and KN channel 
thresholds, to which it couples strongly. The A c (2595) lies below the DN and just slightly above the 7rS c thresholds. 
This similarity was exploited in the first exploratory work [481 ] where free space amplitudes were constructed from 
a set of separable coupled-channel interactions obtained from chirally motivated Lagrangians upon replacing the s 
quark by the c quark. The 1 = A c (2595) was generated as a DN s-wave resonant state of binding energy sa 200 MeV 
with a width of « 3 MeV. While giving the first indication that the A c (2595) could have a dynamical origin, this work 
ignored the strangeness degree of freedom due to its very construction. Therefore, the 7r and K (Goldstone) mesons 
were not treated on an equal footing, and the role of some channels that would appear from the corresponding SU(4) 
meson and baryon multiplets, such as D S A, £> s £j_ KE C and K"E! C , was ignored. A different approach which respected 
the proper symmetries, was attempted in Ref. [49 ] . There, charmed baryon resonances were generated dynamically 
from the scattering of Goldstone bosons off the ground-state J p = i charmed baryons. The C = 1, S = I = 
resonance found at 2650 MeV was identified with the A c (2595) in spite of the fact that, due to the strong coupling 
to the 7rS c states, its width came out to be more than twenty times larger than the experimental width of about 
4 MeV. The trouble with this model is that, due to its very construction, it does not account for the coupling to 
the DN channel which contributes strongly to generating the narrow A c (2595) according to Ref. [48|. A substantial 
improvement came in a recent work [50] in which the above mentioned shortcomings have been overcome by exploiting 
the universal vector meson coupling hypothesis to break the SU(4) symmetry in a convenient and well-defined manner. 
This is done by a t-channel exchange of vector mesons between pseudoscalar mesons and baryons in such a way that 
chiral symmetry is preserved in the light meson sector, while the interaction is still of the WT type. We note that 
in the subsequent sections this approach is referred as TVME (t-channel vector meson exchange model). The model 
generates a narrow C = 1, I = resonance that is identified with the A c (2595), together with other resonances in 
other strangeness-charm sectors with J = 1/2 - . An extension of the model to d-wave J = 3/2~ resonances, generated 
by the interactions of pseudoscalar mesons with baryons of the J = 3/2+ SU(4) multiplet, was developed in Ref. [5l| . 
Finally, some modifications over the model in Ref. [5fJ, both in the kernel and in the regularization scheme, were 
implemented in Ref. [52j, obtaining qualitatively similar results to the earlier approaches. In all these works the 
zero-range limit in the i-channel exchange of vector mesons is identified as the driving force for the s-wave scattering 
of pseudoscalar mesons off the baryon ground states, and is used to carry out a coupled-channel analysis. 

A serious limitation of this SU(4) TVME model is that, whereas the pseudoscalar mesons D and D s are included 
in the coupled-channel dynamics, their vector partners D* and D* are completely left out. This is not just ified from 
the point of view of Heavy Quark Symmetry (HQS), which is a proper QCD spin-flavor symmetry [53, [5J, [55| when 
the quark masses become much larger than the typical confinement scale, Aqcd- 

HQS predicts that all types of spin interactions vanish for infinitely massive quarks: the dynamics is unchanged 
under arbitrary transformations in the spin of the heavy quark (Q) . The spin-dependent interactions are proportional 
to the chromomagnetic moment of the heavy quark, hence are of the order of 1/rriQ. The total angular momentum 
J of the hadron is a conserved quantity, and the spin of the heavy quark Sq is also conserved in the ttlq — > oo limit. 
Consequently, the spin of the light degrees of freedom Si = J — Sq is conserved in that limit. Thus, heavy hadrons 
come in doublets (unless s; = 0, with Sf = si(si + 1)) containing states with total spin j± = Si ± 1/2 obtained by 
combining the spin of the light degrees of freedom with the spin of the heavy quark sq = 1/2. These doublets are 
degenerated in the m,Q — > oo limit. This is the case for the ground state mesons D and D* or D s and D* s which 
are composed of a charm quark with sq = 1/2 and light degrees of freedom with s/ = 1/2, forming a multiplet of 
hadrons with spin and 1 and negative parity. The entire multiplet of degenerate states should be treated in any 
HQS inspired formalism as a single field that transforms linearly under the heavy quark symmetries [54J [55j. For 
finite charm quark mass, the pseudoscalar and vector D meson masses differ in just about one pion mass 3 , even less 
for the strange charmed mesons, thus it is reasonable to expect, for instance, that the coupling DN — D*N might 
play an important role, and could well be bigger than that played by some of the other channel couplings considered 
in the TVME model, say, the DN ~ K c rj. 

The mass splitting between pseudoscalar and vector mesons is comparatively larger in the light sectors, for example, 
between the it and p mesons or the K and the K* . This appears to justify the usual approximation of ignoring vector 
degrees of freedom in the SU(3) coupled-channel studies in s-wave meson (8i)-meson (8i) and meson (8i)-baryon 
(82, IO4) scattering and resonances, within different renormalization schemes. Nevertheless, as already mentioned, 



2 Besides kinematical similarities (open and closed channels near the resonance mass ...), the underlying dynamics might be different as 
a strange quark is related to chiral symmetry and a charm quark to heavy-quark symmetry, as we will discuss below. 

3 One can easily deduce that rn^, — m,£>* = O ( -r ) 



vector mesons have been incorporated in the WT meson-baryon chiral Lagrangian 37, 38, 39, 64] through a scheme 
that treats the six states of a light quark (u, d or s with spin up, j, or down, J.) as equivalent, leading to an SU(6) 
symmetric scheme. Here we will extend this model to SU(8) spin-flavor symmetry to account also for the charm 
quark degree of freedom. Therefore, compared to the TVME approach, our coupled-channel analysis will consider the 
additional role of vector mesons (D*, D*, p, K*,...). Moreover, both models differ on how they implement the strong 
breaking of the SU(4) flavor symmetry observed in nature. 

We should clarify here that the extension of the WT interaction to vector mesons [37[ and to flavor SU(4) is a 
model. Chiral symmetry is a symmetry of the light quark sector, i.e. at best for flavor SU(3), but not for SU(4). 
Chiral symmetry breaking fixes model independently the strength of the lowest order interaction between Goldstone 
bosons and other hadrons (here baryons) - the WT interaction. However, chiral symmetry does not fix the interaction 
between vector mesons and baryons. On the other hand, HQS connects vector and pseudoscalar mesons which contain 
charm quarks. It does not tell anything about mesons made out of light quarks. Nevertheless, it is clearly appealing 
to have a model for four flavors and for pseudoscalar and vector mesons which reduces to the WT interaction in the 
sector where Goldstone bosons are involved and which incorporates heavy-quark symmetry in the sector where charm 
quarks are involved. The model assumption in the present extension does not appear to be easy to justify, but we 
want to try it based on the reasonable outcome of the SU(6) extension in the three-flavor sectors on the one hand, 
and formal plausibility of how the SU(4) WT interaction in the charmed PS meson- SU(3) baryon interaction did 
come out in the vector meson exchange picture as discussed in the TVME approach, on the other hand. We just want 
to fuse the two approaches. Note than in our case, we improve on previous models since we incorporate HQS in the 
charm sector. 

In this exploratory work, we have studied only non-strange single charmed resonances. Our approach should 
work better close to the relevant thresholds, and in particular to describe the lowest lying resonances in each of the 
examined I J (meson-baryon isospin and total angular momentum) sectors. Actually, working near threshold is the 
only justification to ignore the contribution of d-wave meson-baryon interactions and of p- wave couplings, through the 
neglected s- and u-channel diagrams, which otherwise might play an important role. These contributions have been 
also traditionally ignored in the previous successful SU[3) flavor studies of lowest-lying charmless s— and d— wave 
resonances @, 1, 1,1, 0, & S M M 0, 0, M M, M HI M M M M M M,MMMM,MMM, which wc 
aim to extend to the charm sector. 

The present work is organized as follows. In Sec. [II] we gi ye the details of our SU(8) extension of the WT meson- 
baryon Lagrangian. The basis of states and the general form of the Lagrangian are presented in Sec. Ill Al while the 
symmetry breaking effects are discussed in Sec. Ill Bl The unitarization method and renormalization procedure are 
described in Sect. IIIH Our SU(8) results are discussed in Sec. lIVl while Sec.|V]is devoted to a comparison with results 
obtained in a reduced SU(4) model. A summary of our conclusions is presented in Sec. [VlJ 

II. SU(8) EXTENSION OF THE WT MESON-BARYON LAGRANGIAN 

A. SU(8) symmetry 

Since we assume that the pure SU(4) (flavor) transformations commute with the pure SU(2) (spin) transformations, 
it follows that an SU(8) multiplet can be decomposed into SU(4) multiplets each with definite total spin. With 
the inclusion of spin there are 64 quark-antiquark (qq) states, and the irreducible reduction of the SU(8) group 
representation (denoting the SU(4) multiplets of dimensionality n and spin J by n-^j+x) reads 

8 <8> 8* = 63 © 1 = (16i © 15a © la) ©li ■ (1) 

S v ' 

63 

Assuming that the lowest bound state is an s-state, and since the relative parity of a fermion-antifermion 
pair is odd, the SU(4) 15-plet of pseudoscalar (D s ,D,K,Tr,r],r] c ,K,D,D s ) and the 16-plet of vector 
(D*,D*,K*,p, u), J/&, K*,D*,D*, <p) mesons are placed in the 63 representation. Note that the 63 allows nine 
light vector mesons but only eight 0~ light mesons. A ninth 0~ meson must go into the 1 of SU(8). We use pure 
cc wave functions for the charmonium states r\ c and J/^f, and the usual quark content, rj = -4= (uu + dd — 2ss) , 

rf = -4= (uu + dd + ss), u> = -4= (uu + da) , <f> = —ss, for the physical isoscalar light mesons. Such a specification 

induces some mixing between the isoscalar SU(4) mathematical states to build the physical states. Detailed spin-flavor 
wave functions, specifying our conventions are given in Appendix [AJ Mesons of spin J > 1 can be understood as 
states of the qq system with orbital angular momentum L > 0, or molecular-type qqqq meson-meson states. 
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TABLE I: Summary of the quantum numbers and experimental masses of the baryons containing heavy quarks. I, and Si igiLt 
are the isospin, and the spin parity of the light degrees of freedom and 5", J p are strangeness and the spin parity of the baryon 
(I denotes a light quark of flavor u or d). Isospin averaged experimental masses are taken from Ref. 58], with errors counting 
for the mass differences between the members of the same isomultiplet and the uncertainties quoted in Ref. [5q |. 



In the case of baryons, with the inclusion of the spin, one finds 512 three quark states: 

120®56®168®168 = 

(4 4 © 20 2 ) ©2 x (20' 2 © 20 4 © 20 2 © 4 2 ) , 
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where I I and I I I I arc the 20 and 20' SU(4) representations, respectively. It is natural to assign the lowest- 
lying baryons to the 120 of SU(8). This is appropriate because in the light sector it can accommodate an octet of 
spin-1/2 baryons and a decuplet of spin-3/2 baryons which are precisely the SU(3)-spin combinations of the low- 
lying baryon states (N, E, A, S and A, E*, 5*, $1). The remaining states in the 2O2 and 20 4 are completed with the 
charmed baryons: A c , £ c , B c , BJ, ,O c , B cc , fl cc and £*, B*, Q,*, Q,* c , B* c , Q, ccc , respectively. Quantum numbers of the 
charmed baryons are summarized in Table HI The 120 of SU(8) is totally symmetric, which allows the baryon to be 
made of three quarks in s-wave (the color wavefunction being antisymmetric) . 

Here we will focus on the s-wave interaction between the lowest-lying meson (63) and the lowest-lying baryon 
(120) SU(8) multiplets at low energies. Note that at higher energies, higher partial waves would become important, 
so a suitable treatment of spin-orbit effects in the SU(8) scheme should be considered. Now, assuming that the s-wave 
effective meson-baryon Hamiltonian is SU(8) invariant, and since the SU(8) decomposition of the product of the 63 
and 120 representations yields 



63 <g> 120 = 120 © 168 © 2520 © 4752 . 



(3) 



we conclude that there are only four independent Wigner-Eckart irreducible matrix elements (WEIME's) each of which 
being a function of the meson-baryon Mandelstam variable s. The WEIME's might be constrained by demanding 
that the SU(8) amplitudes for the scattering of the Goldstone J p = _ mesons of the pion octet off the J p = | 
baryons of the nucleon octet reduce to those from SU(3) chiral symmetry. At leading order in the chiral expansion, 
these latter amplitudes are obtained from the WT Lagrangian which, besides hadron masses, only depends on the 
pion weak decay constant / ~ 90 MeV in the chiral limit. 

This procedure is a natural extension of the one derived in Ref. [37| for SU(6) to the spin-flavor SU(8) symmetry 
group. If one follows the arguments in this reference, it is clear that not all the SU(3) invariant interactions in the 
(8i)meson-(82)baryon sector can be extended to an SU(8) invariant interaction. This is easily understood because the 
number of independent SU(3) WEIME's to describe the (8i)meson-(82)baryon interaction is six, as follows trivially 
from the SU(3) decomposition 
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whereas the SU(8) spin-flavor symmetry requires the knowledge of four WEIME's, as seen from Eq. ([3]). However, 
Chiral Symmetry (CS) at leading order is much more predictive than SU(3) symmetry, and it predicts the values of 
the various SU(3) WEIME's which otherwise would be totally arbitrary functions of s. Indeed, the WT Lagrangian 
is not just SU(3) symmetric but also chiral (SUj,(3) (g) SUr(3)) invariant. Symbolically, up to an overall constant, the 
WT interaction is 

£ WT = Tr([Mt ! M][Bt )B ]). (5) 

This structure, dictated by CS, is most suitable to analyze in the t-channel. The mesons M fall in the SU(3) 
representation 8 which is also the adjoint representation. The commutator [Aft, Af] indicates a i-channel coupling 
to the 8 a (antisymmetric) representation, thus 

£ WT = ((M t ®M) 8a ®(B t ®B) 8 J 1 . (6) 

Note that from the group point of view there would be as many different t— channel SU(3) singlet Lagrangians as 
WEIME's. 

For the SU(8) spin-flavor symmetry, the mesons M fall in the 63 which is the adjoint representation and the baryons 
B are found in the 120, which is fully symmetric, and the group reductions 

63 © 63 = 1 8 63 s © 63 a © 720 © 945 © 945* © 1232 
120 120* = 1© 63 ©1232 ©13104 (7) 

lead 4 to the total of four different t— channel SU(8) singlet couplings (that can be used to construct s-wave meson- 
baryon interactions) 

((Aft © M) x © (B f © B) 1 ) 1 , ((Af f © Af) 63a © (B f © B) 63 ) 1 , 

((Aft ©Af) 63s ©(B t © B) 63 ) 1 , ((Aft ©Af) 1232 ©(S t © 5)1232)^ (8) 

which match in number to that of the independent WEIME's expected from the group reduction of Eq. (J3]). To ensure 
that the SU(8) amplitudes will reduce to those deduced from the SU(3) WT Lagrangian in the (8i)meson-(82)baryon 
subspace, we set all the couplings in Eq. IjHJ) to be zero except for 

£w4 8) = (( Mt ® M )es a © (S f © B)ea) ± (9) 

which is the natural and unique SU(8) extension of the usual SU(3) WT Lagrangian. 5 The reduction of this Lagrangian 

to the SU(6) sector reproduces the C w j, found in Ref. [371 ]. The corresponding BS approximation, which we will 
discuss below, successfully reproduces the previous SU(3)-flavor WT results for the lowest-lying s- and d-wave negative 
parity baryon resonances obtained from scattering of the mesons of the pion octet off baryons of the nucleon octet 
and delta decuplet [38|. To compute the matrix elements of the SU(6) WT interaction C w ^ , the SU(6)-multiplet 
coupling factors found in [53] were used in [33,(38}. Here, because of the lack of these factors for the SU(8) group, we 
use quark model constructions of hadrons with field theoretical methods to express everything in tensor representations 
as described in Appendix[A] Thus, we get the tree level amplitudes (we use the convention V = —C): 

2 



where the last factor is due to the spinor normalization convention: uu = vv = 1, as in Refs. [1J,[22|]. In the above 
expression IJSC are the meson-baryon isospin, total angular momentum, strangeness and charm quantum numbers, 
M (E) the common mass (CM energy) of the baryons placed in the 120 SU(8) representation, and D IJSC a matrix 
in the coupled channel space. For instance in the J = 3/2, 1 = 2,S = 0,C = 1 sector we consider 5 channels, 
S c/ o, AD, S*7r, AD*, and S*p, without including channels with double charmed baryons since they are much higher 
in energy. The D matrices relevant for this work are presented in Tables IXVIIXXlTII of Appendix [Bj Note that we 
truncate the coupled channel space by not including either double charmed baryon entries, or r\ c , or J/^f mesons, 
since in this work we are just interested in the charm C = 1 sector. Double charmed baryons or ?/ c , or J/^P mesons 
will contribute to this sector only in conjunction to another charmed hadron, and the threshold of the corresponding 
channel would be much higher in energy than those of channels involving a single charmed hadron. 



4 The singlet representation 1 only appears in the reduction of the product of one representation by its complex-conjugate. The couplings 
of Eq. (O arise because 1, 63 and 1232 are self-complex conjugate representations. 

5 Alternatively, the unique spin-flavor symmetric extension of the WT Lagrangian is derived in [64J for arbitrary number of flavors and 
colors assuming only a non-relativistic reduction for the baryons. 



B. SU(8)-flavor symmetry breaking effects 

The SU(8) spin-flavor symmetry is strongly broken in nature. So we have taken into account mass breaking effects 
by adopting the physical hadron masses both in the tree level interaction of Eq. (TIT)]) and in the evaluation of the 
kinematical thresholds of different channels. 

Next, we have considered spin-flavor symmetry breaking effects due to the difference between the weak non-charmed 
and charmed, as well as pseudoscalar and vector meson decay constants. These symmetry breaking effects are sizable 
because, for instance, the ratios fn/fn of fp/fn deviate from one, and actually they are of the order of 1.7. The 
pseudoscalar meson (P) decay constants, fp, are defined by 

(0|gi7 M 7 5 g 2 (0)|P(p)) = -K/2/p^ (11) 

and vector meson (V) decay constants, fv, by 

(0\qi^q 2 {0)\V(p,e)) = VW/v^ , (12) 

where q%, qi are the quark fields, e M is the polarization vector of the meson, and ray its mass. With the above 
definitions, HQS predicts /d (s) = Id* , up to Aqcd/^c corrections 55], which guarantees that the normalizations 
of the coupling constants in Eqs. (TTT]) and (p~2|) are consistent. For light mesons there exist sizable corrections to the 
HQS-type relation fp = fv- 

Within our present approach, besides the use of physical hadron masses, we further break the SU(8) symmetry of 
Eq. (|10p by the replacement: f 2 —> f a fb, depending on the nature of the mesons involved. In Table 1111 we compiled 
the values of the decay constants used throughout this work. Whenever possible, we use the experimental values for 
the decay constants. 

Taking into account all the SU(8) breaking effects stated above, our tree level amplitudes now read 



v usc ( r , P 7J5C 2 V^-AJq-M b j Eg + Mg j E b + M b 

v <* h^ = A* — ifj b — \I^MT\l^mr> (13) 

where M a (M&) and E a (Eb) are, respectively, the mass and the CM energy of the baryon in the a (b) channel. Though 
the way these flavor breaking effects have been introduced is intuitive, we should acknowledge here again that this 
constitutes a further assumption of our model, and that we can not establish a clear connection to any kind of first 
principles, as it is the case for all previous models. 

We finish this section by comparing our present scheme with the one in TVME. The latter relies on SU(4) symmetry 
to construct the effective interaction between pseudoscalar mesons in the 16-plet with the 2O2 and 2O4 plet represen- 
tations of J p = 1/2+ and J p = 3/2+ baryons through a t— channel exchange of the 16-plet of vector mesons. Note 
that vector mesons do not enter in the s-channel coupled space. Then the universal vector meson coupling hypothesis 
provides the global interaction strength among the above S U(4) -multiplets. However, deviations up to 50% from the 
universal value might be expected, as acknowledged in Ref. [50j from an analysis of the D* — > Dtt decay rate. Then, 
aided by the KSFR relation, which is consistent with chiral symmetry at very low energy and momentum transfer, 
the resultant lowest order meson-baryon interaction in the SU(4) limit is found to take the WT form of Eq. (fTOjl in the 
zero range limit [52J 6 . SU(4) symmetry is broken in the TVME model by using the physical hadron masses, without 
adopting different values for meson decay constants. The use of the physical masses for the pseudoscalar mesons and 
the baryons (external legs in the coupled channel formalism of the TVME model) is analogous to the mass breaking 
effects included in the present scheme upon replacing (,/s — M) by (2y/s — M a — M )/2 in Eq. (fTO]) . However, the 
use of physical masses for the vector mesons in the i-channel in the TVME model leads to an additional symmetry 
breaking effect which is different from those induced by the use of different meson decay constants in our scheme. For 
instance, let us consider the E c 7r — > S c 7r, ND — > ND, and E c 7r — > ND transitions. Within the TVME model the first 
two amplitudes are driven by the ^-channel exchange of light vector mesons such as p, while the last one is by the 
exchange of a charmed D* vector meson. Thus, besides different SU(4) coupling coefficients, the last channel obtains 
an additional suppression factor of ~ rri \/m 2 D , = 1/6.8 relative to the first two [52|. Within our scheme, these three 
channels will scale as 1//^, 1//|> and l/(/z)/ir), respectively. Then on top of the matrix elements {D 1 ^ ) we obtain 
suppression factors f^/fp, ~ 1/2.9 and / w //d ~ 1/1.7 for the second and third amplitudes, respectively, relative 
to the first one. In summary, in the present work diagonal transitions involving charmed mesons are suppressed by 



A small p— wave contribution of the TVME model was neglected to reduce it to the WT form [52 



Decay Constant 


Exp. Value [MeV] 


Theor. Estimates [MeV] 


u 


92.4 ±0.1 ±0.3 [58] 


- 


Ik 


113.0 ±1.0 ±0.3 [58] 


- 


fv 


2i 1.2/* [58] 


- 


fp 


~ 153" 


- 


!k* 


~ 153 6 


- 


/«. 


~ 138 c 


- 


n 


~ 163 d 


- 


fo 


157.4 ±11.8 ±2.2 [58, 65] 


- 


f Da 


193.7 ±9.2 ±4.9 [66] 


- 




193.0 ±11.3 ±4.9 [6JJ 


- 




200.1 ±12.0 ±9.9 [68] 


- 


/d» 


- 


165±6±13 [69] c 




- 


223 ±18 [70/ 




- 


164 ±15 [71] 9 


fo*. 


- 


180±4±11 [69] A 




- 


237 ±18 [70]' 




- 


231 ± 13 [71? 


Decay Constant 




Value used in this work 


fv' 




fv 


fn. 




193.7 MeV 


f D , = f D . 




fo 



a From V(p -> e+e~) = 47r» 2 / 2 /3m p and T(r -» p^ T ) = G 2 cos 2 c / 2 (m 2 - m 2 ) 2 (2m 2 + m 2 )/87rm3. 

6 From T(r -> iCV T ) = G 2 sin 2 0c/lr« (m? - m|,) 2 (2m|. » +m 2 )/8vrm3. 

c Prom r(w -» e+e") = 47ro 2 / 2 /27m aJ . 

d From T(0 -> e+e - ) = 87ra 2 /|/27m . 

E Quenched LQCD. 

■'Rolativistic constituent quark model. 

9 Non-relativistic constituent quark model 

h Quenched LQCD. 

'Rolativistic constituent quark model. 

J Non-relativistic constituent quark model 



TABLE II: Meson decay constants. We use V(p -> e + e") = 7.04±0.06 KeV, T(u -> e+e") = 0.60±0.02 KeV, T(. 
1.27 ± 0.04 KeV, T(r -> pz/ T ) = (5.71 ± 0.07) x 10~ 10 MeV and T(r -> K*v T ) = (3.06 ± 0.45) x 10" 11 MeV from 
G = 1.1664 x 10" 11 MeV -2 and a = 1/137.036 are the Fermi and fine-structure constants, and cos#c = 0.974, 
the Cabibbo angle. 



i-»eV) = 
581 ] . Besides, 
the cosine of 



about a factor of 3 with respect to the TVME SU(4) models [50, [54], while the charm-exchange amplitudes will not 
be suppressed by the mass factor ~ m 2 D , /m 2 = 6.8 but only by a factor of about two due to the use of physical values 
for the decay constants. Thus, it is clear that the difference in the pattern of SU(4) symmetry breaking in the two 
schemes discussed above is even qualitative. 

III. BS MESON-BARYON SCATTERING MATRIX 

We solve the coupled channel BS equation with the interaction kernel determined by Eq. (|T5|) which incorporates 
the SU(8) breaking effects discussed in Subsect. Ill Bt For a given IJSC sector, the solution for the coupled channel 
s-wave scattering amplitude T IJSC (^fs) 7 in the on-shell scheme @, H S 122, ISa, ISa, IfiO] is. 



^usc 



(VS) = 



1-V IJSC (^)G IJSC (^) 



V 



use 



(VS). 



(14) 



Here G IJSC (y/s) is a diagonal matrix consisting of loop functions. The loop function for channel i reads, 

Ol'» WS) = i2A, / |L {p _ q)2 [ M? + m -^^ , (15) 

where Mi and m; are the masses of the baryon and meson in the channel, respectively. This quantity is logarithmi- 
cally divergent, hence must be regularized. For example, the regularization may be done by one-subtraction at the 
subtraction point y/s = n{ JSC such that 

G^ c (v^ = M?' SC ) = 0, (16) 

with index i running in the coupled channel space. The TVME model of Refs. [50, uLU] uses the Renormalization 
Scheme (RS) of Eq. (|16p . and we adopt here its choice for the value of the subtraction point /j,l JSC . It is taken to be 
independent of the spin J, then set identical within a given sector ISC to \Jrtx^ + M^, where m t h + A/th is the mass 
of the lightest hadronic channel. For the SU(6) sector, this RS recovers previous results for the lowest-lying 1/2 - and 
3/2~ baryon resonances appearing in the scattering of the octet of Goldstone bosons off the lowest baryon octet and 
decuplet [3, @, ffl , flfj. Ha . [23 . [23l . [23, y0|, and leads to new predictions for higher energy resonances. According to the 
authors of Refs. [50L J5l|. such a choice guarantees an approximate crossing symmetry although such a claim appears 
somewhat dubious because crossing symmetry involves isospin mixtures, thus choosing an alternative subtraction point 
might lead to yet another reasonable result. Nevertheless, one should bear in mind an apparent correlation between 
the renormalization procedure and/or subtraction point, and the values for the meson decay constants used in the 
potential V IJ s . For instance, it should be useful to mention the SU(6) results which are somewhat complementary 
to the present approach in the choice of the subtraction points and meson decay constants Ref. |38| | . By taking a 
different value of the subtraction point for each channel i of a given IJSC sector in the form \i\ JSC — \J ' Mf + rnj , 
the model gives a good agreement with a common Goldstone boson decay constant value of 90 MeV (pion decay 
constant in the chiral limit), instead of using the different values quoted in Table HT1 

Given these uncertainties, in the present work we allow for slight changes in the choice of the subtraction point by 
introducing a parameter a, such that 

(// sc ) 2 = a(m t 2 h + M t 2 h ) , (17) 

which will be fitted to data. In principle, a could have a different value in each IJSC sector. Unfortunately we 
have insufficient empirical information. Therefore, we will adjust the value of a to reproduce the position of the well 
established A c (2595) resonance with IJSC = (0, 1/2,0, 1), then the same value will be used in all other sectors. 

The mass and widths of the dynamically generated resonances in each IJSC sector are determined from the positions 
of the poles, zr, in the second Riemann sheet of the corresponding scattering amplitudes, namely Mr = Re {zr) and 
Tr = — 21m (zr). The coupling constants of each resonance to the various baryon-meson states are obtained from the 
residues by fitting the amplitudes to the expression 

Tl^iz) = g '. e g ' e , (18) 

(z - zr) 

for complex energy values z close to the pole, where the complex couplings are written in terms of the absolute 
value, <7fc, and phase, <fik- These dimensionless couplings determine the corresponding decaying branching ratios to the 
open baryon-meson channels. Since the dynamically generated states may couple differently to their baryon-meson 
components, we will examine the jj-channcl independent quantity 

f^%)^ mfx J2\T^%)\, (19) 

i 

which allows us to identify all the resonances within a given sector at once. 

IV. SU(8) RESULTS 

In this work we have focused on the non-strange (S = 0), singly charmed (C = 1) baryon resonances. The 
experimental status in this sector is summarized in Table UTTl In the following tables and figures, we collect results of 
resonances up to 3.5 GeV, which is basically the energy of the heaviest channels used in the present approach. Note, 
however, that the position and width of dynamical states generated ~ 500 MeV beyond the lowest threshold for each 
IJSC sector must be taken with caution. In reality, they will be strongly influenced by other ingredients not included 
in the present model, such as the coupling to three-body meson-meson-baryon states or interactions of higher angular 
momentum. This observation applies especially to those states higher than 3 GeV that couple strongly to low lying 
baryon-meson channels. 



Resonance I(J P ) Status Mass [MeV] F [MeV] 



A c (2595) 


0(1/2") 


*** 


2595.4 ± 0.6 


3.6 ±2.0 -1.3 


A c (2625) 


0(3/2") 


SfCSfclfC 


2628.1 ±0.6 


< 1.9 


A c (2765) 


?(? ? ) 


* 


2766.6 ± 2.4 


50 


or E c (2765) 










A c (2880) 


0(5/2+) 


SfCSjCJjC 


2881.9 ±0.5 


5.8 ±1.9 


A c (2940) 


0(? ? ) 


*** 


2939.8 ± 1.6 


18 ±8 


E c (2800) ++ 


1(? ? ) 


*** 


2801 ±4-6 


75 ± 22 - 17 


£ c (2800)+ 


1(? ? ) 


*** 


2792 ±14-5 


62 ± 60 - 40 


E c (2800)° 


1(? ? ) 


*** 


2802 ±4-7 


61 ± 28 - 18 



TABLE III: Summary of experimental data for baryon resonances with charm 1, zero strangeness and negative (or unknown) 
parity as compiled in Ref. [53]. The A c and E c ground states are omitted. 



0, J = 1/2 



In this sector the following 16 channels are involved: 



£ C 7T ND A c r] ND* E C K A c u E' C K AD S 

2591.6 2806.15 2833.91 2947.27 2965.12 3069.03 3072.52 3084.18 

AD* S E cP A cV ' E*p A c E C K* E' C K* E* C K* 

3227.98 3229.05 3244.24 3293.46 3305.92 3363.33 3470.73 3540.23 



where the second line gives the channel thresholds in MeV. 
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FIG. 1: The function T = maxj J± \Tij\ in the I = 0, J — 1/2 sector calculated along the scattering line within our broken 
SU(8) model, for different values of the subtraction point. The stars denote the nominal position of the A c (2595) resonance. 



In Fig. [T] we show our results for T as a function of the CM. energy y/s, where we have used the experimentally 
known meson decay constants and have chosen fp a — 193.7 MeV. For rj', D* and D*, we have taken f,y = f v and 
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1=0, J=1/2, a =0.9698, f„=f„ =f n 

s.p. D D D 

s 




2400 2500 2600 2700 2800 2900 3000 3100 3200 3300 3400 

Re s 1/2 [MeV] 



FIG. 2: The function f = max., £). \Tij\ in the I — 0, J — 1/2 sector, calculated within the present broken SU(8) model 
with a = 0.9698. Stars show the positions of the resonances obtained, while open cycles denote those of the experimentally 
measured resonances, the A c (2595) in this figure. 

TABLE IV: Mass [MeV], width [MeV] and couplings of resonances with S = and C = 1, generated in the / = 0, J = 1/2 
sector using our broken SU(8) model. See text for more details. The third column quotes only the values of the most important 
couplings, and those shown in bold correspond to channels open for decay. 



/ = 0, J = 1/2 


Mr 


T R 


Couplings to main channels 


2595.4 


0.58 


9Ec „ = 0.36, g ND = 3.69, g ND * = 5.70, g ADs = 1.42, g AD% = 2.94 


2610.0 


70.9 


g Sc7r = 2.25, g ND = 1.47, g ND * = 1.81, g ScP = 1.22 


2821.5 


1.0 


9jvd = 0.32, g AcV = 1.2, g~ cK = 1.79, g AD ; = 1.11, gs cP = 1-23, g s * p = 1.15 


2871.2 


91.2 


gND = 2.0, g A D s = 1.15, g N D> = 2.15, gAD; - 1-92, g AcU = 1-01, gs cP = 2.56, g S ; P = 0.94 


2937.2 


71.7 


9A c r7 = 1.34, g A D 3 = 1.4, g N D* — 1.51, c/ad* = 3.41, g s * p = 2.23 


2954.7 


65.4 


gs^ = 1.02, g^ K = 1.2, g A o s = 0.85, g A ^ = 2.46, g^ p = 1.16, gs iP = 0.91 


3020.1 


47.0 


gs c K = 1.13, g ADs = 1.07, jadj = 1.46, g^ cP = 1.51, g s * P = 2.49, g~ tK * = 0.95 


3053.3 


3.49 


9A cU = 0.11, 9 ad s = 1.43, g s , rK = 1.49, g Ac <t, = 1-27, ffs'p = 1.02 


3152.1 


94.4 


9ad„ = 1.74, ffADj = 2.02, g Ac4 , = 2.17, g KK * = 0.98, 5Sc a- = 1.38, <?h;k* = 1.16 


3244.7 


87.4 


9ad„ = 0.72, g AD * = 0.74, g SaK = 0.68, gs / K - = 2.32, g^ K » = 2.63 


3272.6 


4.31 


9s; . K = 0.17, g ScP = 0.17, g A D„ = 0.15, g Ac ^ = 1.37, g Bc K* = 2.26 


3315.9 


43.4 


Ss cP = 0.47, g A D B = 0.62, g Ac * = 0.66, g Hc K* = 1.2, ffa'„x* = 1.91, 3h;a- = 2.38 



/rj» = /d* = Jd (see Table HI)) . The dotted line shows that a narrow resonance is produced very close to the position 
of the nominal A c (2595) which is denoted by the starred line in the figure. By slightly changing the value of the 
subtraction point in the renormalization scheme, viz. setting the scaling factor a to 0.9698 [see Eq. (|17p], we can 
reproduce the position of the A c (2595). 

There are several other narrow resonances and bumps in Fig. [T] which can be better disentangled by investigating 
the behavior of the function T and T in the complex plane. We obtain the poles of T in the second Riemann sheet 
and determine the coupling constants to the various baryon-meson channels through the residues of the corresponding 
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TABLE V: Positions of poles (in MeV) generated in the J = 1/2, 1 = sector from uncoupled calculations in the present 
broken SU(8) model with a = 0.9698. An asterisk * denotes resonances. 



Channel (M thrcs ) 


2600 


2700 


2800 


2900 


3000 


3100 


3200 


3300 


3400 


£ c tt (2592) 


2621.6* 


















ND (2806) 




2775.4 
















Ac-q (2834) 




















ND" (2947) 


2687.3 


















E C K (2965) 








2947.7 












A c lo (3069) 




















E' C K (3073) 










3042.8 










AD S (3084) 




















AD* (3228) 










3057.3 










E c p (3229) 








2962.6 












A c r]' (3244) 




















S*p (3293) 








2988.5 












A c <f> (3306) 




















H CJ ft"* (3363) 














3293.3 






3^* (3471) 
















3382.6 




E* C K* (3540) 


















3440.2 



amplitudes, as indicated in Eq. (fl~8"|) . The positions of the poles of T are visualized in Fig. [5] as peaks of T, while 
the values of the mass and width of the corresponding resonances, together with the values of the couplings to the 
baryon- meson components to which they couple most strongly, are collected in Table IIV1 

The width of the A c (2595) resonance turns out to be small: 0.58 MeV, smaller than the experimental width, 



3.6 



-2.0 
-1.3- 



However, we have not included here the three-body decay channel A c 7T7t which already represents almost 



one third of the decay events [58| . The narrowness of this resonance is due to the fact that it lies only a few MeV 
above its unique two-body decaying channel S c 7r to which it actually couples very weakly, as indicated by the low 
value of the corresponding coupling constant in Table ITVl 

Here a remarkable feature is that there is a second resonance very close to A c (2595). This is precisely the same 
pattern found in the charmless I — 0, S — —1 sector for the A(1405) |22| . 

It is illustrative to compare the coupled channel results with those from an uncoupled calculation as shown in 
Table. |Vj In the lower energy region we find poles at 2622 MeV, 2687 MeV and 2775 MeV, which are, respectively, 
bound S c 7r, ND* and ND states. Note that the bound states in ND and ND* channels are quite close. This 
emphasizes once more the importance of including the vector mesons in the description. The effect of the channel 
coupling may be easily traced by taking into account the largest coefficients in Table IXV1 together with the channel 
threshold energies. We conclude that the resonance at 2595.4 MeV appears to be basically the remnant of the ND* 
bound state found at 2687 MeV in the uncoupled calculation, but with the important additional binding effects from 
the ND channel, to which it couples very strongly, along with some moderate modifications induced by coupling to 
AD S and AD* the latter lying more than 400 MeV away. As emphasized in the next section, we observe here that 
the nature of this resonance is substantially different from that found in Refs. [5fJ, [52j , where it is identified as a ND 
bound state. A wider resonance at 2610 MeV originates from the S c 7r resonance found at 2622 MeV in the uncoupled 
calculation, mildly modified by its coupling to ND as well as to some channels involving vector-mesons such as ND* 
and Yj c p. We do not find a clear remnant of the uncoupled ND state at 2775 MeV, apart from the already commented 
influence on the properties of the A c (2595) resonance. 

A narrow resonance at 2822 MeV originates mainly from the E^if bound state at 2947 MeV, which is substantially 
influenced by the coupling to the A c r; channel as well as to a few other channels (AD*, E cj o and S*p) involving 
vector mesons. The two other narrow resonances found in this sector, at 3053 MeV and 3273 MeV, can be traced 
back, respectively, to the "E! C K bound state at 3043 MeV and the S C K * bound state at 3293 MeV of the uncoupled 
calculation. The resonance at 3053 MeV shows also a strong coupling to a few other channels, especially to the 
neighboring AD S state. 

Most of the remaining wider resonances show also a stronger coupling to a baryon-meson system for which the 
uncoupled calculation produces a bound state. They are simply wider because they also couple significantly to 
meson-baryon states lying below their mass. However, we identify two resonances that occur genuinely as a result of 
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the coupled channel formalism, since their main baryon-meson component (A c lu for the resonance at 2955 MeV and 
A c (j) for the resonance at 3152 MeV) does not have a bound state in the uncoupled calculation. 

In order to study possible variations of the result due to input parameters, we have also performed calculations by 
adopting the values for the unknown decay constants of the charmed vector mesons from a quenched lattice QCD, 
LQCD (69J: f D . = 165 MeV, f D * = 180 MeV. Then the A c (2595) is reproduced with a slightly smaller scaling factor 
(a = 0.9558), but the qualitative features for the positions, widths and couplings of all the other resonances generated 
in this sector do not change drastically. 

In comparing with data, we look at the possibility of identifying some of our narrow states with a resonance seen 
experimentally. Our predictions are limited by the fact that we have implemented neither the coupling to three- 
body states nor p-, d- or higher-multipolarity interactions. In the present I = 0, J = 1/2 sector, we have already 
identified our first narrow state with the A c (2595) and have established it as being a quasi-bound ND* system. A 
resonance at 2880 MeV with a width of T ~ 6 MeV, decaying into A c tttt states with some fraction of resonant decay 



through E c 7r states, was reported in [42J. One might identify it with our narrow resonance at 2822 MeV, especially 
because if this state were moved to the experimental energy of 2880 MeV, it would appear above the threshold of the 
A c rj channel to which it couples significantly. However, the Belle collaboration determined recently the spin of the 
A c (2880) to be J = 5/2 from the E c (2455)7r decay angular distribution, and the parity positive from agreement of 
the r(E c (2520)7r)/r(£ c (2455)7r) branching ratio with a prediction of HQS [zl [7J uA, UM ■ These assignments do not 
match the spin-parity 1/2 - of our dynamically generated state. 

There is yet another resonance in this sector, the A c (2940) of width T ~ 18 ±8 MeV which has been seen in the Dp 
invariant mass distributions, but its spin and parity have not been determined [4y|. If it were a J = 1/2 - state, one 
might be tempted to identify it with one of our narrow states. However, none of the resonances obtained by our model 
having a width smaller than 50 MeV, except for the lowest-lying one already identified with the A c (2595), couples 
significantly to the ND states from which the invariant mass of the A c (2940) is reconstructed. Our model does not 
give either any resonance in this energy region that couples preferentially to ND* states, the threshold of which is 
only a few MeV above the A c (2940) mass. Therefore, we do not expect the A c (2940) to be a molecular ND* system, 
as claimed in the literature [76[ . We finally note that the pD° histogram shown in Ref . [46] is not incompatible with 
the existence of a very narrow state just above the ND threshold MeV, as the one at 2822 MeV and width 1 MeV 
found in the present work, as well as in the dynamical model of Ref. [501 ] . 

B. 1 = 1, J = 1/2 

In the I = 1, J = 1/2 sector there are 22 channels involved: 

A c tt S c tt ND ND* E C K Z c rj A c p E' C K Y>D S AL>* E c p 

2424.5 2591.6 2806.15 2947.27 2965.12 3001.01 3061.95 3072.52 3161.65 3218.35 3229.05 

S c w £*p S*w T.D* E C K* Y, c r)' E' C K* E c £*£>* £*0 S*AT* 

3236.13 3293.46 3300.54 3305.45 3363.33 3411.34 3470.73 3473.02 3496.87 3537.43 3540.23 

Our results are presented in Fig. [3] and Table [VT] First, we find a very narrow resonance at 2554 MeV, which is 
basically composed of the AD* component (somewhat mixed with the ND and T,*D* states). Hence, this resonance 
would be absent in models not including channels with a vector meson and a 3/2 + baryon. Kincmatically it can only 
decay to A c 7r, but the corresponding coupling is very weak. 

At higher energies and up to 3.4 GeV, our model generates many resonances in this sector, most of them having 
rather appreciable widths. Some of these resonances find their analogue within the states obtained by the pioneer 
work of Lutz and Kolomeitsev |49j , where only the meson-baryon channels composed by an octet Goldstone boson 
plus a member of either the anti-triplet or the sextet of open charm ground state baryons were considered. The 
resonances at 2612 MeV, 2637 MeV and 2974 MeV, would correspond to the states in Ref. [49|] found, respectively, at 
2800 MeV (coupling mostly to A c tt), 2700 MeV (coupling mostly to £ c tt) and 2985 MeV (coupling mostly to E c rj). 
The states obtained here appear at a lower mass and show a different width due to their coupling to many of the new 
baryon-meson components considered in this work. 

From the experimental side, an isotriplet of excited charmed baryons, £ c (2800) has been reported [431 ]. decaying 
mainly to A+7r~, A+7r° and A+7r + pairs with a large width of around 60 MeV (but with a > 50% error). This 
resonance has been tentatively assumed to decay to A c 7r in d-wave with its spin parity assumed as J p = 3/2~ relying 
on quark models, but angular distributions have not yet been measured. Of the predicted states discussed above 
having a moderate width, only the one at 3096 MeV shows a sizable coupling to A c ir states, but it lies 300 MeV 
away of the experimental mass. Then what about identifying the predicted 3035 MeV state as the E c (2800)? Note 
first, that this 3035 MeV state decays significantly to S c 7r pairs. Now if this state were to move to the experimental 
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1=1, J=1/2, a =0.9698, f=f=f n 
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FIG. 3: The same as Fig. [2]for the I = 1, J = 1/2 sector. 



position, it would become somewhat narrower, especially because two of the possible decay channels, E c t? and ND* , 
would be closed. Its decay to A c tt pairs, as observed experimentally, would require an implementation in the present 
model of an additional meson-baryon interaction in <i-wave or/and p-wave coupling in the neglected s- and w-channel 
diagrams. However, according to the sizable coupling to S c 7r states, the identification of our state as £ c (2800) would 
be ruled out if the later were not observed in experiments looking at the A c 7T7t systems. On inspecting the histograms 
of Ref. [Hj], one finds signals around M(A+7r+7r~) - M(A+) ~ 500 MeV, corresponding to the mass of the £ c (2800), 
although they are probably too feeble. 

It seems unlikely that any of the three states with energies lying between 2800 and 3000 MeV could qualify for the 
observed isotriplet, basically because they couple quite significantly to ND states. Hence, if the parameters of our 
model would allow us to move one of these states below the ND threshold, it would become significantly narrower 
than the experimental width of ~ 60 MeV. However, on inspecting the histogram of I = 1 D + p pairs, shown in 
Ref. [46| to confirm the A c (2940) as being an isosinglet, one observes a clear enhancement around 2860 MeV of width 
~ 10 MeV, which we could identify with one of the states we see decaying into ND pairs. The most likely candidate 
is the state at 2974 MeV because, if it was to be moved to 2860 MeV, its width would be reduced from 37 MeV to a 
value of around 10 MeV due to the closing of two out of three decaying channels. 

Above 3 GeV we find a few states that c oup le strongly to some vector meson-baryon channels. We thus expect 
them to be absent in the SU(4) models [5CJ, [52]. It may be of interest to note that the relatively narrow resonances 
at 3158 and 3346 MeV could be seen from the invariant mass of T, c r] and £D S states, respectively. 



2, J = 1/2 



In the I = 2, J = 1/2 sector there are 4 channels involved: 



£ C 7T AD* S c/ 9 E*p 

2591.6 3218.35 3229.05 3293.46 



As seen from the results shown in Table IVTIl our model generates two narrow structures, at 3172 MeV and 3246 
MeV which qualify, respectively, as E c p and S*p quasi-bound states. 
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TABLE VI: The same as Table [TVl for the J = 1, J = 1/2 sector. 



7 = 1, J = 1/2 


Mr 


r fl 


Couplings to main channels 


2553.6 


0.67 


g AcW = 0.15, g ND = 2.28, g AD * = 6.74, g s * D , = 2.89 


2612.2 


179.0 


9A a T! = 1-95, g ND . = 3.78, g EcP = 1.27, g E ; p = 1.4 


2637.1 


79.9 


9s C 77 = 1-98, g ND = 2.35, 5jvz>* = 1.69, # e *b; = 1-24 


2822.8 


34.8 


9jvd = 1-55, g SDs = 1.01, g Se K - 1.04, grjvr>* = 1.41, #a cP = 1-37, g A D* = 1.81,g E > = 2.27, g-j/t* = 


0.94 


2868.0 


38.6 


gND = 0.76, 0jvx>* = 2.75, g ED . = 1.05, g AcP = 1.3, p EcP = 1-54, g SaUJ - 1.02, p A D* = 1-78 


2974.0 


37.2 


g ND = 0.72, g Sc7r = 0.63, gjVD . = 0.63, g ScV = 1.67, <? AcP = 1.29, g ScLJ = 1.13, <?ad* = 1.21 


3001.9 


73.8 


<?s c *r = 1.42, g SDs = 1.68, ffED: = 4.15, 5EcP = 1.42, <? E . p = 1.26 


3034.9 


64.2 


gs c „ = 0.72, g Scr? = 0.67, g^D* = 0.76, g^ Ds = 1.16, g Sc u = 2.32, ffE . D j = 2.44, g s> = 1.16, ffs; 


_, =0.87 


3077.7 


48.7 


gjVD . = 0.85, g SclJ = 2.40, gs-o* = 3.20, ff S * w = 1.43 


3095.9 


44.0 


<?a c „ = 0.39, g s , nK = 0.88, gjv ,D* = 0.45, g SDs = 1.34, 5EcP = 1.09, ff S » DJ = 2.80, gv*u> = 2.36 


3147.5 


48.8 


<M cP = 1-02, 3EcP = 1.62, ffScU = 0.94, g s> = 1.67, g s> = 1.57 


3158.2 


8.6 


9s C 77 = 0.3, g AcP = 0.27, g EDa = 1.54, # ED * = 1.51, g~ a K* = 1-53 


3265.8 


36.8 


9sd„ = 0.69, g SD * = 1.54, p Ec0 = 1.32, g^ K , = 1.49, gx^ = 1.67, 3h -k- = 1.34 


3316.2 


20.6 


SS-.K = 0.44, ffEc0 = 1.66, gape. = 2.22, flE . z ,. = 1.80, 5h; k* = 1.01 


3346.1 


9.0 


9sd 8 = 0.45, ffEc0 = 1.81, gs c K- = 1.12, ffE ^ = 1.38, 5 h-k- = 1.14 


3392.8 


15.5 


g s . p = 0.31, <? E . = 2.17, g s ?K* = 1.95 
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FIG. 4: The same as Fig. [2] for the 7 = 2, J = 1/2 sector. 



D. 



7 = 0, J = 3/2 



Of various J = 3/2 resonance candidates, we first consider those in the 1 = sector, with 11 coupled channels: 



S*tt ND* K c u E*K AD* E c p E*p A c </> E C K* E' C K* E* C K* 
2656.01 2947.27 3069.03 3142.02 3227.98 3229.05 3293.46 3305.92 3363.33 3470.73 3540.23 
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TABLE VII: The same as Table |TVl but for the / = 2, J = 1/2 sector. 



I = 2, J = 1/2 


Mr 


Tr 


Couplings to all channels 


2936.5 


185.0 


9s C 77 = 1-62, g AD * - 5.37, g ScP = 


1.12, g s> = 0.66 


3171.5 


8.26 


9Se77 = 0.33, g AD , = 1.61, g EcP = 


2.26, ffE . p = 0.76 


3246.1 


8.09 


9s C 77 = 0.18, g A D* = 0.33, g ScP 


= 0.47, g^, p = 2.39 
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FIG. 5: The same as Fig. [2] for the I 
resonances, A c (2625) in this case. 



0, J = 3/2 sector. Open cycles denote the positions of the experimentally measured 



The results displayed in Fig. [5] and Table IVITlI show that, in this sector, our model generates nine resonances, four 
of which have a width smaller than 5 MeV. 

Experimentally, there is a three star A c (2625) resonance with J p = 3/2~, which has a very narrow width, T < 1.9 
MeV, and decays mostly to A c 7T7t [58|]. It has been identified as the counterpart of the A(1520) which, in dynamical 
models, shows a strong coupling to E*7r [29|, [30, [3l|, [38|, [30] ■ In the present S = 0, C = 1 sector, we indeed find a 
I = 0, J = 3/2 d-wave A c resonance at 2660 MeV with its width r = 38 MeV, which couples very strongly to the S*7r 
channel. We may therefore identify this resonance as a strong candidate for the d-wave A c (2625). A small change 
in the subtraction point fi IJSC could easily move the resonance down by 40 MeV to the nominal position. By so 
doing it should also be possible to reduce the width considerably, hence closer to the experimental value because the 
resonance position would get below the threshold of the channel E*7r to which this resonance couples strongly. We 
note that a similar resonance was found at 2660 MeV in the SU(4) model of Ref. [5l\. The novel feature in the recent 
approach is that the resonance in question here has non-negligible baryon-vector meson components such as ND* 
and Tj* c p. Other narrow resonances observed in this channel are a state at 3094 MeV coupling strongly to Z*K, as 
in Ref. [5l| , but also to the vector- meson channel AD* and three narrow resonances which couple mostly to channels 
with vector meson components. Out of the resonances obtained in this sector around 3 GeV and beyond, the one 
at 2941 MeV could be a candidate for the A c (2940) of width ~ 18 MeV and unknown J p , observed recently from 
the invariant mass of D°p pairs 46] . Its decay into D°p pairs would imply the additional implementation of p-wave 
interactions in our model. The strong coupling of this resonance to ND* hints also for a substantial coupling into 
p-wave ND states. 
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TABLE VIII: The same as Table ITVl for the I = 0, J = 3/2 sector. 



I = o, J = 3/2 


Afe 


Th 


Couplings to main channels 


2659.5 


37.8 


9s «„ = 2.23, g ND . = 2.11, g s> = 1.34 


2940.5 


2.06 


g s . w = 0.21, g ND . = 2.21, g s> = 1.40 


2972.8 


64.9 


g s , w = 1.05, ffAcW = 2.42, ffE - K = 0.95, ffE . p = 1.28 


3093.5 


3.48 


9A cU = 0.29, g ND * = 0.28, <? ao . = 2.89, g^ K = 1.88 


3180.5 


13.4 


ffB . K = 0.69, g AD ; = 2.49, g Ac4 , = 1-83, flSoK * = 1-00, g B * K * = 


0.79 




3209.8 


0.6 


9JVD . = 0.08, g B . K = 0.10, 9EcP = 1.81, gv>„ = 1-09 


3274.1 


3.48 


9ad; = 0.19, 9B ;k = 0.2, g Ac ^ = 1.30, <?s c x* = 2.31 


3330.4 


36.9 


g Ar>i = 0.50, g Aa <t, = 0.68, g s . p = 0.57, g ScK " = 0.85, ffs^K* 


= 1.85, g Si K* 


= 2.34 


3396.3 


4.8 


9ad s . = 0.17, g Ac = 0.22, g s . p = 0.15, 9s / K , = 2.17, <?h;k* 


= 1.86 





E. 1 = 1, J = 3/2 

In the 1 = 1, J = 3/2 sector one finds 20 channels: 

£*7T ND* A cP £*?/ AL> E*K AD* E c p E c w E*p 

2656.01 2947.27 3061.95 3065.42 3077.23 3142.02 3218.35 3229.05 3236.13 3293.46 

T,*u> £L>* £*L> S S CJ RT* S^if* S c £*r/ £*£>* £*</> S*i<f* 

3300.54 3305.45 3353.07 3363.33 3470.73 3473.02 3475.75 3496.87 3537.43 3540.23 



100 
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FIG. 6: The same as Fig. [2]for the I = 1, J = 3/2 sector. 

Our results are presented in Fig. [5] and Table ITXl We obtain many resonances and they all couple strongly, without 
exception, to channels with vector mesons. A bound state at 2550 MeV, whose main baryon-meson components 
contain a charmed meson, lies below the threshold of any possible decay channel. It would develop a narrow decay 
width to A c 7r states if we incorporated d-wave interactions. This resonance is the charm sector counterpart of the 
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TABLE IX: The same as Table ITVl but for the / = 1, J = 3/2 sector. 



I = 1, J = 3/2 


Mr 


r R 


Couplings to main channels 


2549.8 


0.0 


g ND * = 2.52, g SD . = 2.22, g AD = 4.23, g s * Ds = 1.48, g AD , = 5.28, g s * Di = 2.29 


2686.9 


60.0 


g s .„ = 1.91, g ND * = 2.68, g SDt = 0.96, g AcP = 1.02, g^» D , = 0.97 


2901.6 


14.7 


g s .^ = 0.58, g ND , = 2.77, g AcP = 1.52, ffA D = 1-26, g AD > = 1.05, ffE . p = 1.20 


3012.1 


22.7 


g ND , = 0.90, g SD * = 2.05, <? E ^ = 2.04 


3026.5 


11.2 


g ND . = 0.72, ffEc „ = 2.10, g AD = 2.40, g AD . = 2.83 


3062.6 


25.5 


9s^77 = 0.61, g-£ D * — 1.67, g A D* = 1-56, g E *o; = 1-69, p E « w = 1.97 


3125.1 


4.05 


g^ = 0.17, fls .„ = 0.21, ffEO; = 1.39, g^ Da = 1.67, g s * K = 1.33, 5e . D j = 2.90, 


<? E ^ = 0.93 


3156.1 


36.1 


9A cP = 0.9, gED j = 1.45, g ScP = 1.50, g ScU = 1.27, g E *p = 1.41, gx* u = 1.39 


3179.5 


8.89 


9s^ = 0.37, g-znt = 2.65, g SaP = 1.40, ffHc K« = 1-71, </s> = 0.85 


3196.4 


17.3 


9a cP = 0.56, gsD' - 0.96, gz aP = 1.17, gs c ^ = 0.87, g~ c K' = 0.87, g S j P = 2.16, g S ' 


w = 0.92 


3277.2 


0.62 


gs * K = 0.11, g Sc4 , = 0.91, te / x . = 1.68, <7e-d 3 = 2.61, flE » D j = 2.52 


3345.4 


19.1 


gs , K = 0.54, g siK * = 1.18, flE » D « = 1.07, g s ^ = 1.38, ga. K . = 2.57 


3352.6 


6.65 


9ed; = 0.40, g s * K = 0.32, <? Ec0 = 2.26, g^ K , = 1.10, g ScK * = 0.91, <; E ^ = 1.72, 


g s , c K* = 0.99 


3386.3 


7.79 


g ScK , = 0.48, ffEc0 = 1.40, g KK * = 1.00, ffE ^ = 2.20, ffS jK- = 1.18 



hyperonic S(1670) resonance, which has also been seen in chiral unitary models [29l. I30L l3l| coupling strongly to the 
AK channel. Indeed, the charmed bound state found here couples strongly to the corresponding AD state in the 
charm sector, and even more strongly to its vector partner: AD*. This is precisely the reason for the large binding 
of this resonance, compared to that of its strange sector counterpart. 

Since the observed 1=1 charmed baryon resonances do not have spin nor parity assignments, most of the arguments 
discussed in the / = 1, J = 1/2 section, can also be made here. We might assign the resonance at 2687 MeV, decaying 
primarily to S*7r, to the bump at 2760 MeV observed from A c irir states |42j. However, the corresponding width 
estimated by scaling our result: T2760 = (<7276o/<?269o)r2690, would turn out to be around 120 MeV, much larger than 
the experimental value. Then, similarly to our discussion in the I = 1, J = 1/2 sector, the next state at 2901 MeV, 
decaying into E*7r s-wave pairs, could be a candidate for the E c (2800), seen in A c tt pairs [43], if this resonance were 
also seen in A c Tm states. The little enhancement in the A c tttt histogram [42| around M(A+7r + 7r~) — M(A+) ~ 500 
MeV, which is obviously distorted by the excitation of the lower mass resonance at 2760 MeV, might well correspond 
roughly to the mass of the S c (2800). 

Our model predicts a few narrow resonances above 3 GeV. Of particular interest are those coupling strongly to 
ND* , since they could be observed in DnN invariant mass distributions and, if sufficiently narrow, these states should 
be reflected into the DN invariant mass spectra Although the discussion in Ref. [46| mentions explicitly the absence 
of signal in the D*(2010)+p or D*(2007)°p mass distributions, the invariant mass distribution of D + p pairs shows 
some structures which, unfortunately, are embedded into the statistical error fluctuations. Clearly, more data needs 
to be collected in order to see whether some of these structures acquire statistical significance to become new charmed 
isotriplet baryons. 



F. 1 = 2, J = 3/2 
Here we present our results in the I = 2, J = 3/2 sector, where there are only five channels: 



£*7T AD AD* S cj o S*p 

2656.01 3077.23 3218.35 3229.05 3293.46 

Our results, presented in Fig. [7| and Table IXl indicate that the model generates a state at 2911 MeV, with mainly 
ND and ND* components, which is very wide because it lies above the threshold of S*7r states to which it also couples 
significantly. We also find a narrow resonance at 3071 MeV, which lies only a few MeV below the AD threshold and 
coupling more strongly to AD* states, and a wider resonance at 3196 MeV which qualifies as a T,*p quasi-bound 
state. 



18 



1=2, J=3/2, a =0.9698, f =f n , =f_ 

s.p. D D D 

s 




2850 2900 2950 3000 3050 3100 3150 3200 

Re s 1/2 [MeV] 



3250 



FIG. 7: The same as Fig. [2] for the I = 2, J = 3/2 sector. 



TABLE X: The same as Table |P7l but for the I = 2, J = 3/2 sector. 



J = 2, J = 3/2 


Afe 


r R 


Couplings to all channels 


2910.8 


158.0 


gs*7i- = 1.67, qad — 3.52, gAD* — 3.95, gs cP — 


0.72, ffs * p = 


0.68 


3070.8 


9.79 


Se'tt = 0.40, g A D - 1-48, pad* = 2.73, gs oP = 


1.35, gs* P = 


1.18 


3195.8 


25.7 


gs*7i- = 0.50, g/AD = 0.51, gAD* = 1-07, gs cP 


= 0.97, flfc.. 


= 2.38 



G. J = 5/2 (J = 0, 1 and 2) 



We finally present our results for the isospin sectors 7 = 0, 1 and 2 with maximum value of spin J = 5/2. The 
states building up the coupled-channel space in this J = 5/2 case are composed of a 3/2 + baryon and a vector meson. 
Note that, as seen from the coefficients in Tables IXXIl IXXLTl and IXXIIIl all these cases present the peculiarity that 
channels without strange quarks (as S*p or AD*) do not mix with channels having an associated ss pair (as £*</> or 
S*AT*). Therefore, in the case I = 0, J = 5/2, since the interaction is strongly attractive in the two channels involved, 
namely 

Kp z* c k* 

3293.46 3540.23 , 

we obtain two states at 3258 MeV and 3440 MeV that are, respectively, bound S*p and E*K* systems, as seen in 
Table |XH Since we do not take into account three-particle channels, these resonances appear as stable bound states. 
In the 7 = 1, J = 5/2 sector there are 6 channels involved: 

AD* S*p S*w Z*D* £*</> Z*K* 
3218.35 3293.46 3300.54 3496.87 3537.43 3540.23 



Table IXIII shows that we generate 3 bound states and a narrow resonance. The lower two states couple to the 
non-strange baryon-meson channels, while the upper two belong to the associated strangeness sector. 
Finally, we present in Table IXlTTl the results for the I = 2, J = 5/2 sector, which involves two channels: 
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TABLE XI: The same as Table ITVl but for the I = 0, J = 5/2 sector. 



I = o, J = 5/2 


Mfl 


r fl 


Couplings to main channels 


3258.0 


0.0 


g s * p = 2.24 


3440.2 


0.0 


g^ K , = 2.87 



TABLE XII: The same as Table \N\ but for the / = 1, J = 5/2 sector. 



1=1, J = 5/2 



Mr 



3084.5 



3233.8 



3276.7 



3421.9 



0.0 



0.17 



0.0 



0.0 



Couplings to main channels 



gAP- = 3.99, g E ; p = 0.98, ffs;^ = 1.92 



g Ag » = 0.14, ff S » p = 2.14, ffs^ = 3.14 



ffs»p* = 4.32, g H » = 1.17, g s * K * = 1-95 
= 0.15, gs;0 = 2.52, g^ K « =1.65 



gs-o; 



AD* £*p 
3218.35 3293.46 



TABLE XIII: The same as Table ITVl but for the I = 2, J = 5/2 sector. 



I = 2, J = 5/2 


M B 


r fl 


Couplings to all channels 


3126.6 


0.0 


gAD' = 3.62, g s * p = 2.0 



As seen from the coefficients in Table IXXIIIl the interaction is zero for AD* states and repulsive for S*p states. In 
spite of this, our model generates a bound state at 3127 MeV, about 90 MeV below the threshold of the AD* channel 
and coupling strongly to it. This bound state is entirely due to the coupling between AD* and S*p states, facilitated 
by the significant value of the mixing coefficient and by the vicinity of the channels involved, whose thresholds are 
only 75 MeV apart. We have checked that the bound state disappears if the mixing coefficient is reduced to 2/3 its 
actual value. 



V. COMPARISON WITH SU(4) MODELS 



In this section we would like to compare our results with those obtained within a standard SU(4) model denoted as 
TVME in the present work. We will follow the TVME approach, first developed in Ref. [50(, but with modifications 
implemented in Ref. [52} , which we describe briefly in the following. Firstly, the last term of the interaction kernel 
used in Ref. [50] , namely 



vv> s > c HVS)=J2 



c 



(I,S,C) 
V 



8 ml 

ye [16] v 



(2VI- 



M -M+ (M - M) ■ 



(20) 



is dropped. This was shown in [52| to be more consistent with the zero range (t — > 0) limit of the TVME employed in 
Ref. [50l | . In the second place, a common heavy vector meson mass ttih and a common light vector meson mass ray 
is used in the kernel, with the relation ran ~ 2mv ■ In this way, the resulting interaction is of the same form as that 
shown in Eq. (|13p . namely 



V. 



use 

ab 



(vG)=XcDi, 



4/2 



usc 2 Vs-M a ^M b jE a + M a jE b + M, 



2M a 



2M h 



(21) 



but with an additional factor \c, which is \ c = 1/4 for charm-exchange transitions and x c = 1 otherwise. This factor 
accounts approximately for the ratio (my /tub) 2 between the squared masses of uncharmed and charmed vector 
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mesons, such as the p and D* , respectively. The additional scalar-isoscalar attractive interaction characterized by a 
^dn term, implemented in Ref. [52j , is not considered here, since it would not affect the qualitative aspects of our 
results. 

The main difference between the SU(4) models [49|, [50, lU, [52| and our SU(8) model, which is a must to be consistent 
with heavy-quark symmetry, is that the latter includes the vector meson-baryon channels. Another essential difference 
lies in the transition amplitudes between states containing heavy mesons because, according to Eq. (|13j) . they scale 
with the inverse of a heavy-meson decay constant for each heavy meson involved, whereas in the SU(4) models the 
decay constant is kept fixed for all transitions. Therefore, as discussed at the end of Sect. Ill Bt amplitudes such as 
ND — + ND or AD S — > ND are typically a factor 3 smaller than those of the SU(4) models. 

TABLE XIV: The same as Table ITVl but for the SU(4) model, see more details in the text. 



Mr 


r R 


Couplings to all SU(4) channels 


1 = 0, J =1/2 


2595.4 


2.01 


g So7r = 0.67, g ND = 6.03, g AcV = 0.12, gs c K = 0.07, g S ' r K = 0.17, g ADs = 3.08, ff Ac y = 0.29 


2625.4 


103.0 


9Se7r = 2.30, g ND = 1.55, g AaV = 0.04, g EaK = 0.03, g B i K - 0.67, g ADs = 1.05, g AcV > = 0.1 


2799.5 


0.0 


g So7r = 0.35 ■ 10" 2 , g ND = 0.05, g AaV = 1.47, g~ aK = 2.57, ffS '„K = 0.02, g A D s = 0.26, g AaV < = 0.02 


3024.8 


31.3 


9s e ^ = 0.59, g ND = 0.50, g Acrl = 0.16, g ScK = 0.11, g^ K = 2.22, g ADs = 1.48, g AcV , = 0.02 


I =1, J = 1/2 


2661.2 


36.5 


g A ^ = 0.65, g So7r = 1.1, g ND = 4.87, gs c „ = 0.04, g ScK = 0.48, g B ' r K = 0.45, g^ Ds = 3.99, gv cV > = 


0.35 


2694.7 


153.0 


9A C 7r = 0.34, g So7r = 2.35, g ND = 1.95, gs cV - 0.21, g~ cK = 0.31, g 3 > r K = 1.19, gsD s - 2.03, g ScV > = 


= 0.18 


2938.0 


291.0 


g A ^ = 1.87, g Sc7r = 0.09, g ND = 0.32, j Ec , = 0.13, g Sc K = 2.46, ffB ^ K = 0.07, g SDs = 1.35, (? Ec ^ 


= 0.19 


2956.5 


29.5 


9A C 77 = 0.05, gs e ^ = 0.63, gND = 0.55, gs C 7, = 1.85, gs c K = 0.26, g s > r K — 1-60, gED a = 1.55, g SaV 


= 0.03 


J = o, J = 3/2 


2678.0 


82.7 


fls> = 2.23, ffH .K = 0.61 


3118.7 


29.6 


gs » w = 0.67, g B * K = 2.06 


I = 1, J = 3/2 


2543.5 


0.0 


Sb.« = 0.22, g^ = 0.12, g AD = 7.21, (? s; k = 0.07, 3 e*d s = 3.27, ff E « v = 0.28 


2745.5 


155.0 


g s;7r = 2.2, 5E;?7 = 0.22, pad = 0.45, g sf K = 1.09, 3 e*d s = 0.52, g sw = 0.06 


3042.4 


31.2 


gs,, = 0.74, ffE .^ = 1.75, g AD = 0.04, 3H jk = 1.52, g s * Ds = 0.07, ff s «^/ = 0.02 



The results obtained within the SU(4) approach, Eq. (|2"Tj) , are collected in Table IXTVl only for those IJ combinations 
which generate resonances. That excludes the 7 = 2, J = 1/2 sector, which has a repulsive interaction in the only 
allowed channel, S c 7r; as well as the two-channel (S*7r, AD) 1 = 2, J = 3/2 sector, where the mixing coefficient is 
not strong enough to overcome the repulsive E*7r interaction. Note that we have slightly modified the subtraction 
point with a parameter a = 0.959, see Eq. (jTTJ), in order to place the A(2595) resonance on the right energy. The 
generation of the resonance would have been facilitated if an attractive scalar-isoscalar Hdn term were included, in 
which case the subtraction point would have to be modified by a value of a = 0.979, closer to one. 

The main observation in the I = 0, J = 1/2 sector is that the model is able to generate a narrow resonance 
around 2600 MeV, which can be readily identified with the A c (2595). It couples strongly to ND and its nature is 
therefore very different from the one obtained in the SU(8) model, which qualifies as being mainly a ND* quasibound 
state. The reason for the SU(4) model being able to generate the A c (2595) is that the lack of baryon-meson states 
in this sector (namely, those involving a vector meson and/or a decuplet baryon) is somehow compensated by the 
enhanced strength in the baryon-meson channels involving a charmed meson. In other words, a resonance that has a 
strong ND* component in the SU(8) approach may reappear as a ND state in the reduced SU(4) model space due 
to the enhanced interaction kernel. We also find relatively narrow resonances at 2799 and 3025 MeV that couple, 
respectively, mostly to S c if and "BI C K states involving uncharmed meson channels and, therefore, are the analogous 
states as those obtained within the SU(8) model at 2821 and 3053 MeV respectively [see Table HV] . although the last 
state in the SU(4) model couples appreciably to S c 7r and ND, and hence is substantially wider. Logically the SU(4) 
model does not produce the resonances that couple strongly to states with light vector mesons in the SU(8) approach. 

In the I = 1, J = 1/2 sector, we observe that the SU(4) model obtains a narrow resonance at 2661 MeV which 
couples strongly to ND. Since the strength of this amplitude is reduced, the SU(8) model does not find such a type 
of resonance, but instead, a very narrow one at 2554 MeV that couples strongly to AD* and more moderately to 
£*£)* and ND, see Table |VH So again, while the two resonances do not differ that much in energy, their very nature 
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is substantially different. The other relatively narrow resonance at 2957 MeV, coupling mainly to E c 7y and Z' C K, 
finds its analogue with the one at 2974 MeV obtained in the SU(8) model, although this one has also vector-meson 
components of comparable size. The narrow states at 3158 MeV and 3346 MeV, which couple mainly to vector meson 
components and decay into E c ?7 and £D S pairs, respectively, are not obtained in the SU(4) approach. 

In the / = 0, J = 3/2 sector, the SU(4) model produces resonances at 2678 MeV and 3119 MeV that couple 
mainly to £*7r and S*iif, respectively. The first one is equivalent to that found at 2659 MeV in the SU(8) model, see 
Table. IVIIII The second one could be identified basically with that appearing at 3094 MeV in the SU(8) model, that 
also has stronger vector-meson components. It has been argued in previous works [49 . |5l| and in the present one, 
that the lowest-lying resonance could correspond to the observed three star A c (2625) of J p = 3/2~, which is quoted 
to be the analogous in the charm sector of the A(1520) [58j. The SU(8) model produces, apart from the state at 2940 
MeV that could be identified with the A c (2940), a few other narrow resonances beyond 3 GeV that couple strongly 
to states containing vector mesons and, therefore, cannot be seen in the SU(4) model. 

In the I = 1, J = 3/2 sector, there are three resonances that couple mainly to AD, £*7r and S*?y, respectively, and 
are the analogies of some of the lower energy states obtained in the SU(8) model, namely those at 2550 MeV, 2687 
MeV and 3012 MeV. Note, however, that the resonance at 2550 MeV, which is the counterpart in the charm sector 
of the hyperonic £(1670) resonance having mainly AK components [23 . |30 . l3lj . becomes a AD* resonance in the 
SU(8) model, although its coupling to AD, ND*, £D*, and £*D* is also significant. The three resonances obtained 
here with the simulated SU(4) model are the same as those found in Ref. [51(, however with a notable exception for 
the width of the state around 3 GeV which is quoted to be very narrow, while we find it to be around 30 MeV, due 
basically to a quite sizable coupling to £*7r states. Note that, in this sector, the WT coefficient is —2 for the £*7r 
transition and for the T,*r), so the sizable coupling of this resonance to this latter channel is necessarily due to the 
non-perturbative coupled channels processes involving non-diagonal transitions. 

We end this section by noting that the present simulation of the SU(4) model of Ref. [49] reproduces all the 
resonances given in Refs. [49, (51j at approximately the same positions. Only the widths of some resonances appear to 
be wider here. This is because we are using an interaction kernel that ignores the tensor term of the t-channel vector 
meson propagator, more consistent with the limit t — > 0, which is implicit in the contact WT interaction [52| . 

VI. CONCLUSIONS 

In the present work, we have studied charmed baryon resonances within a coupled channels unitary approach that 
implements, for the first time, the characteristic features of heavy-quark symmetry, as for instance the fact that D and 
D* mesons have to be treated on an equal footing. This is accomplished by extending the t-channel vector-exchange 
SU(4) models (TVME) used in Refs. [49, ISO, [5l|, [52J to SU(8) spin-flavor symmetry, then by implementing a somewhat 
different way of breaking the flavor symmetry. More concretely our tree-level s-wave WT amplitudes are obtained 
not only by adopting the physical hadron masses, but also by introducing the physical weak decay constants of the 
mesons involved in the transitions. This procedure reduces considerably the strength of diagonal amplitudes involving 
a charmed meson, as comparing to SU(4) models. 

The present SU(8) model generates dynamically the resonances with negative parity in all the isospin-spin sectors 
that one can form from a s-wave interaction between the mesons of the _ , 1~ multiplets with the l/2 + , 3/2 + baryons. 
In this work we have focused only on the strangeness and charm quantum numbers appropriate for the ND interaction, 
namely S = and C = 1. In the I = 0, J = 1/2 sector, the lowest narrow state obtained can be readily identified 
with the A c (2595) which in our approach may be interpreted as being primarily a ND* bound state. We may also 
identify our state at 2660 MeV in the I = 0, J = 3/2 sector, which couples strongly to the £*7r channel, to be the 
three star A c (2625). The latter is usually referred to as being the charm counterpart of the A(1520), which shows 
a strong coupling to £*7r in dynamical models. Similarly, in the I = 1, J = 3/2 sector we find a resonance at 2550 
MeV that couples strongly to AD and AD* states. This appears to be the counterpart of the £(1670) hyperon which 
couples to AK in dynamical models. We also find in this sector a state at 2901 MeV that could be a candidate for 
the £ c (2800). 

Our results have been compared with those obtained in the SU(4) models [43, [50, [5l|, [52|]. It appears that our 
SU(8) approach reproduces all the resonances generated in the SU(4) models which couple strongly to the channels 
consisting of a pseudoscalar octet meson and a charmed baryon. They may be identified as chirally excited charmed 
baryons. On the other hand, due to the different pattern of flavor symmetry breaking, resonances in the SU(4) model 
that couple strongly to baryon-meson states containing a charmed meson and an un-charmed baryon (such as ND and 
YD S ) are, in principle, not necessarily reproduced within our SU(8) approach which has adopted the physical decay 
constants, resulting in a reduced interaction in the corresponding transition amplitudes. However, the additional 
states implemented by the enlarged model space in SU(8) due to the requirement of heavy-quark symmetry, such 
as ND* and AD*, compensate largely for the reduced interaction, producing resonances in the same region or even 
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more bound, than those obtained in the SU(4) model, but with a quite different composition. A prime example is 
the A c (2595) in the 1 = 0, J = 1/2 sector. Within the SU(4) models this resonance is dynamically generated mainly 
from ND states and it is then interpreted as a ND quasi-bound state, with a mixture of AD S components, which lies 
very close to the threshold of its only decaying channel, S c 7r. Instead, the SU(8) model interprets this resonance as 
being mainly a ND* quasi-bound state, although with sizable ND, AD S and AD* components. 

Our SU(8) approach predicts more states than those obtained within the SU(4) model. Not all of them find a direct 
identification with one of the observed resonances in the C = 1, S = sector. We note, however, that many of our 
states couple weakly to the baryon-meson pairs from which the resonances are measured. It should also be pointed 
that some of our states would become wider or even disappear in a more realistic model containing also three-body 
channels and higher multipolarity interactions. We finally note that the experimental spectra show a limited amount 
of counts and, in order to disentangle new resonant structures from data, more statistics is clearly needed. 
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APPENDIX A: HADRON SPIN-FLAVOR WAVE FUNCTIONS AND TENSOR METHOD 

This Appendix is aimed for a detailed account of how to calculate the matrix element D^ in the expression for 
the WT interaction, Eq. (13). It consists of defining the hadronic wave functions from the quark model and the tensor 
representation of the SU(8) operators and wave functons. 

1. Wavefunctions for the ground state hadrons 

Pseudoscalar mesons: 
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Spin 1/2 baryons: 



\f>) 



\K* 
\K* 
\D* 

\d; 

\D* 
\D* 



m = 



oi T oi T io), 

-QtiQ%\o)> 

oi T oi T io>, 

-Ql T Qt T |0), 
-Qt r Qlt\°)> 

■±=(qKqU 

-Ql T Qt T |o>, 



Q$M T ) lo> , 



|A> = 

\N) = 

IS) = 

IS) = 

|2c> = 

K) = 

m = 71 

|Ac) - 71 

|Scc> = 71 



1 

71 
i 

71 
i 

71 
i 

7s 

i 

7s 

i 

Ve 

i 

i 

71 

i 

72 

1 

7s 

i 

73 



Qu]Qd[Qs] ~ QuiQd\Q s y I") > 

< y«T < *«T < *<i| _ ^u-\QuiQd]J |0) ' 

Qu^Qs^Qsi - QuiQs-\Qsi 



(A2) 



Q U ]QsiQc] ~ QuiQs]Qc] ) I™ ' 

( ^«T^cX < ycT — ' ^mI^cT^cI,) 1^/ ' 



(A3) 



24 



Spin 3/2 baryons: 
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[Notes] 



• Each state represents a member of a spin and isospin multiplet. The wavefunctions made explicit are those with 
the highest value of 73 and J 3 in the multiplet. 

• Q u -t creates a u quark with spin up, Q \. creates a quark d with spin down, and so on. As usual, color takes care 
of the fermionic statistics, so the quark operators are bosonic. The basis {| f), I I)} la of the standard SU(2) by 
which step operators J± have non-negative matrix elements. Also standard are {\u}, \d)}, {\d), \s}}, {\s}, |c)}, 
and {— \d), \u)}, {— \s), \d)}, {— |c), \s)}. 7 This means that the flavor-SU(4) step operators act as 

E 12 \u) = \d), E 23 \d) = \a), E 34 \s) = |c), 
E 12 \d) = -\u), E 23 \s) = -\d), E u \c) = -|g). (A5) 

• No special choice of relative phases has been made between different SU(4) multiplets. The relative phases 
of the hadronic states in an SU(4) multiplet are taken such that £12, E 23 and E34 have non-negative matrix 
elements. Exceptions are the neutral mesons 77, 77', rj c , ip, 10 and <p. 

• The pseudoscalar meson r\ is (up to a sign) the isoscalar meson of the SU(3) octet and 77' the SU(3) singlet. r\ c 
is purely cc. In the spin 1 sector, ip is purely cc whereas (j) and w have been rotated in the SU(3) sector so that 
<j) is purely ss (ideal mixing). 

2. Tensor representation for matrix elements 

Having identified the group structure of the extended WT interaction, the calculation of matrix elements can be 
done using standard group theoretical techniques. Methods based on computing SU(n) Clebsch-Gordan coefficients 
(e.g. .57]), extracting isoscalar factors (in SU(3) 63], and then SU(3)-scalar factors in SU(4), etc) are useful for an 
explicit hand calculations for small groups but become more involved as the groups get larger. Instead, we follow a 
different route here. A tensor representation spanned by quark operators is used. This method only requires to know 



7 This is the convention used in |6ll, 162(1 for arbitrary SU(n). Note that the convention in ]63| for SU(3) is, instead, {\u), \d)}, {\u), \s)}, 

{-\d), |o>}, {-|g>, |«». 
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Kronecker deltas and Wick contractions. Although it may not be practical for hand calculations, it is conceptually 
simple so that it is easy enough to automatize the calculation for larger groups. 
Following J64J the WT Hamiltonian with spin-flavor symmetry takes the form 

- £wt = ftwT = - ^ : [*, d^ 3 B\ lm B^ m : (A6) 

where B*-* is the baryon field, $ 8 j the meson field, and the SU(8) way of labeling the spin-flavor indices goes for 
i, j, etc. from f to 8 for four flavors. B lJ is completely symmetric under permutation of indices. There is no 
universally accepted convention for the normalization. Here we choose that B l ^ k is normalized in a standard manner 
for a fermionic field when the indices i, j, k are all different. Upon extracting the kinematical factors {y/s— M )/(2/ 2 ), 
the spin- flavor dependence is all contained in the matrices D^. These are the matrix elements of the operator (see 

M) 



{M\M^ k - M^ k M k J )BUe m B^ m : (A7) 

ijk 



where B l i k and Br... denote the baryon annihilation and creation operators, normalized as 



l BiJk > B lj'k'} = S i> S j> S k> + 444 + • ' ' (3! permutations) . (A8) 

Likewise M % j and AiT J 'j = {M % jy denote the meson annihilation and creation operators, normalized as 

[M i j ,M^ j 'ff] = S i i ,S 3 r . (A9) 

These characterizations automatically guarantee that the normalization and spin-flavor transformation properties 
of the states (0\B ljk and (0\M z j coincide with those of (0\Q 1 Q j Q k and (0\Q l Qj, respectively. Here Q l and Qi are the 
SU(8) labeled annihilation operators of quarks and antiquarks, respectively. They transform as 

Q l -* V % 3 Q\ Qi -► (U^YQj = QjiU-^i, U e SU(8) (A10) 

and have standard normalization 

[Q\Q i j] = 5 i j , [Qi,QU] = S>. (Afl) 

The correspondence between the SU(8) and the explicit spin-SU(4) flavor representations are the following: since 
Qi transforms as Q\ , the conventions in subsection IA f I imply the identifications (up to a global sign) 

Q l = {Qu],QdhQshQchQuiiQdi,Qsi,Qci) , 

Qi = {Qui,QdiTQsiiQci,-Qu], -Qdp—Qsf,~~Qci) ■ (Af2) 

That is Q 1 — Q u f, Qs = — Qdp Q^ 6 = — Qj t , etc. The correct book keeping between an explicit flavor-spin index and 
the SU(8) one (running from 1 to 8) needs to be maintained to apply the Kronecker deltas in (|A8|) and (|A9j) . Using 
these identifications together with 

q\q)q\ - 4fc> 

Qt*Q} -> Aft*,-, (A13) 

the mesonic and baryonic wavefunctions can be rewritten in terms of the operators M^j and BJ- k . For instance, 

\N) = -L(Bl 16 -Bl 52 )\0), (A14) 

(once again for the highest weights in the spin-isospin multiplet). ft is straightforward now to obtain the matrix 
elements (MB\Q\M'B r ) using the commutation relations, or equivalently Wick's theorem. A non-trivial confirmation 
of the calculation is to see that the matrices D 1 ^ so obtained have the correct eigenvalues 6, 0, —2, — f 6, —22 (the 
being the consequence of the inclusion of the SU(8) singlet among the mesons). Another important point to check is 
that all the matrix elements computed are (up to a sign) square roots of rational numbers. While this is not surprising 
for matrix elements of SU(n) with standard bases, it is found to be the case also for the neutral vector mesons with 
ideal mixing. 
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APPENDIX B: COEFFICIENTS OF THE s-WAVE TREE LEVEL AMPLITUDES 

This Appendix gives the coefficients D^ sc of the s-wave tree level Baryon-Meson 8 amplitudes of Eq. (fTU|) for the 
various IJSC sectors, in the particular case of S = and C = 1. 



We say Baryon-Meson amplitudes to stress that the given £) IJSC matrix elements corresponds to the state of the two colliding hadrons 
coupled in this order: Baryon-Meson. Change of phases in the matrix elements can be derived of using other order. 
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TABLE XVI: I — 1, J — 1/2, S — 0, C — 1. Baryon-Meson states with more than one c quark have not been included. 
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-2 








f~4 

V 3 


fs 
V 3 

















fs 

V 9 


fi 

V 3 





2 
3 





2 
3 





v» 


-*#* 


v/2 


/T 








-V^ 


V2- 


V 3 





pF 

V 3 





2 
3 


pi 

V 9 


/l28 

V 9 


8 
3 


2 
3 








fs 

V 9 


2 
3 


/s" 

V 9 


/s" 

V 9 


10 

3 



to 

-x 



TABLE XVII: / = 2, J = 1/2, S — 0, C — 1. Baryon- Meson states with more than one c quark have not been included. 



E C 7T 

AD* 

S*p 



E C 7T 

2 

V8 



AL>* 




£ c p 



S*p 



to 



TABLE XVIII: I — 0, J — 3/2, S — 0, C — 1. Baryon-Meson states with more than one c quark have not been included. 



E> 
iVD* 

A c cu 

E c p 
S*P 

A c 

s c ir 
Kk* 

E*K* 



E*7T 

-4 

-V2 



-y/3 







V3 

-1 

-V5 



ND* 

-V2 







A c w 


Ve 



1 



2 

^20 



-1 



-V3 



1 

-2 

-v/2 

-1 

-V5 

V2 



4 
3 

20 
3 



AL> S * 







-V2 







-2 

V2 



E c p 



2 
-1 



_8 
3 

20 
9 



1 



S*p 



/20 
-V5 


20 
9 

_ 16 

3 





A c 




-2 







-y/2 



E C K* 


-1 



1 
-y/2 



S'tf* 



-y/5 





o 



TABLE XIX: I — 1, J — 3/2, S — 0, C — 1. Baryon-Meson states with more than one c quark have not been included. 





E*7T 


ND* 


A c p 


S*?7 


AD 


Hc# 


AD* 


E c p 


E c o; 


E*p 


E*w 


ED* 


E*D S 


s c /r 


H^* 


E c (/) 


s*7 


e*d: 


E*^ 


E*K* 


e> 


-2 


— /i 


-78 





1 


-^ 


/ ^ 

V 3 


[a 

V 3 





pm 
V 3 











72 


— V 3 














pm 
V 3 


ND* 


fi 

V 3 


4 
3 


7? 


fi 
V 9 


— /iS 





/io" 

V 9 


2 
3 


— V 9 


pm 
V 9 


pm 
V 9 


4 
3 


pi 

V 3 











2 
3 


pm 

V 9 








A c p 


-78 


76 











-1 








[a 

V 3 





pm 
V 3 








1 


pi 
V 3 














pi 
V 3 


Kv 





pi 

V 9 








pi 

V 3 


-V3 


/To" 
V 9 














pa 
V 9 


— V 3 


71 


-1 








pm 
V 9 





-75 


AD 


1 


pa 

V 3 





fi 

V 3 


-4 





/so" 

V 3 


pi 


— /s 


V 3 


pm 
V 3 


pm 
V 3 


-2 











pi 
V 3 


_ /is 








Z* C K 


-72 





-1 


-73 











pi 

V 3 


_ pi 


pm 
V 3 


/ ^ 
V 3 


m 
V 3 


v^ 


2 





pi 
V 3 





/is 


— /i2 





AD* 


V 3 


/so" 
V 9 





pa 
V 9 


/io" 
V 3 





20 
3 


pm 
V 9 


pm 
V 9 


1 
3 


pi 
V 9 


/so" 

V 9 


pm 
V 3 











pm 
V 9 


10 

3 








E c p 


V 3 


2 
3 








— /J 


— /i" 


/20" 

V 9 


2 
3 


pa 
V 9 


pm 

V 9 


pm 
V 9 








pi 

V 3 


pi 

V 9 














pm 
V 9 


E C W 





pi 

V 9 


fi 

V 3 





— /i" 


/T 

V 3 


/To" 
V 9 


[32 
V 9 


4 
3 


pm 
V 9 


pm 
V 9 








pi 

V 3 


1 

3 














/ ^ 
V 9 


E*p 


/Is 


pM 

V 9 








V 3 


/To" 
V 3 


1 
3 


pm 
V 9 


pa 
V 9 


10 

3 


pm 
V 9 








/Is 


pm 
V 9 














pm 
V 9 


E*u, 





pa 
V 9 


pX) 

V 3 





/To" 
V 3 


pi 

V 3 


vf 


pm 
V 9 


pm 
V 9 


pm 
V 9 


4 
3 








V 3 


— V 9 














5 
3 


ed: 





4 
3 





/i" 

V 9 


/Is 


[2 

V 3 


/so" 

V 9 














4 
3 


pi 
V 3 


v'fi 


pi 
V 9 


2 
3 


2 
3 


pm 
V 9 


pm 
V 9 


pm 
V 9 


E*D S 





fi 

V 3 





— /i" 


-2 


v^ 


_ /20" 














pi 
V 3 


-1 





p% 
V 3 


pi 

V 3 


pi 


/ ^ 


pi 
V 3 


pm 
V 3 


h c a:* 


72 





1 


V3 





2 





\/f 


pi 
V 3 


pm 
V 3 


\/l 


V6 











pi 

V 3 








— /iS 





KK* 


pi 

V 3 





/T 

V 3 


-1 











pi 

V 9 


1 

3 


pm 
V 9 


V 9 


pi 
V 9 


V 3 





4 
3 


pa 
V 9 





pm 
V 9 


pm 
V 9 


pm 
V 9 


E c (f> 

















pi 

V 3 

















2 
3 


pi 

V 3 


pi 
V 3 


pm 
V 9 








pm 
V 9 





pm 
V 9 


Kv' 





2 

3 








/T" 





pm 
V 9 














2 
3 


pi 














7 








E * D * 





p20 

V 9 


• 


/To" 
V 9 


/20" 
V 3 


pa 
V 3 


10 

3 














pm 
V 9 


~ V 3 





pm 
V 9 


pm 
V 9 


/ ^ 
V 9 


5 
3 


1 
3 


fi 
V 9 


E*0 

















/To" 

V 3 

















pm 
V 9 


-a/| 


pm 
V 3 


pm 
V 9 








1 
3 





fi 
V 9 


-*.#* 


— /i2 





~~ V 3 


-V5 











pa 
V 9 


— V 9 


pm 
V 9 


5 
3 


pm 
V 9 


/iS 





pm 
V 9 


pm 
V 9 





pi 
V 9 


— V 9 


4 
3 



TABLE XX: I — 2, J — 3/2, S — 0, C — 1. Baryon-Meson states with more than one c quark have not been included. 



E*7T 

AD 
AD* 

E c p 



E*7T 

2 



AD 

Vs 





-2 

v/5 



AD* 






E c/ 9 



10 
3 



^P 

V 3 



to 



TABLE XXI: 7 = 0, J — 5/2, S — 0, C — 1. Baryon-Meson states with more than one c quark have not been included. 

S^p -2 

■*K* -2 



co 



TABLE XXII: I — 1, J — 5/2, S — 0, C — 1. Baryon- Meson states with more than one c quark have not been included. 

AD* E*p S*cj E*D* E*0 H*^* 

AD* 2^/8 

E*p 2 Vs 

E*oj \/8 V8 2 

S*DJ -2 \/8 

E*0 0-2 -a/8 

H*AT* \/8 -\/8 2 






TABLE XXIII: 1 = 2, J — 5/2, S = 0, C = 1. Baryon-Meson states with more than one c quark have not been included. 

AD* S*p 

AD* \/l2 

E*p 712 4 



in 
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